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Mechanosensitivity of the 11“^^ kind: 

TGF-/3 mechanism of cell sensing the substrate stiffness 
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Cells can sense forces applied to them, but also the stiffness of their environment. These are two different 
phenomena, and here we investigate the mechanosensitivity of the IP^ kind: how the cell can measure an 
elastic modulus at a single point of adhesion - and how the cell can receive and interpret the chemical signal 
released from the sensor. Our model uses the example of large latent complex of TGF-/d as a sensor. Stochastic 
theory gives the rate of breaking of latent complex, which initiates the signaling feedback loop after the active 
TGF-;d release and leads to a change of cell phenotype driven by the a-smooth muscle actin. We investigate 
the dynamic and steady-state behaviours of the model, comparing them with experiments. In particular, we 
analyse the timescale of approach to the steady state, the stability of the non-linear dynamical system, and 
how the steady-state concentrations of the key markers vary depending on the elasticity of the substrate. We 
discover a crossover region for values of substrate elasticity closely corresponding to that of the fibroblast to 
myofibroblast transition. We suggest that the cell could actively vary the parameters of its dynamic feedback 
loop to ‘choose’ the position of the transition region and the range of substrate elasticity that it can detect. 
In this way, the theory offers the unifying mechanism for a variety of phenomena, such as the myofibroblast 
conversion in fibrosis of wounds and lungs and smooth muscle cell dysfunction in cardiac disease. 


I. INTRODUCTION 

How cells sense is important to life: homeostasis ne¬ 
cessitates sensors. The relationship between cell mor¬ 
phology and chemical properties of its environment have 
long been understood and documented. A plethora of 
sensory systems embedded in the cell membrane are able 
to detect specific signaling molecules, giving information 
to the cell to which it can react. Over the previous two 
decades there has been growing evidence to suggest that 
cells use these chemical pathways to sense the mechanical 
properties of their environment. 

Mechanosensing is a process where at least one chem¬ 
ical reaction in the cell changes in response to a change 
in its mechanical environment. We specifically wish to 
distinguish two very different cases: when external forces 
are applied to the cell, and when the cell needs to mea¬ 
sure the stiffness of its passive environment. We call 
the former case the ‘mechanosensitivity of the kind’. 
There is a lot of research on this problem: the cell re¬ 
sponse to hydrostatic pressure^substrate topography^ 
are just some examples; practically, there are implica¬ 
tions for biotechnology^ and cancer research^’^. Ulti¬ 
mately this kind of mechanosensitivity is well under¬ 
stood: there are many molecular processes that respond 
to an applied force and several different ones are em¬ 
ployed by cells in different situations. It has been ob¬ 
served (for example) that mechanotransductive pathways 
can link compressive forces to translational and post- 
translational events by cytoskeletal interactions^, that 
tyrosine kinases are part of the pathway of the cell re¬ 
sponse to shear stress^, and that mechanical forces can 
alter ion channel permeability^. In each of these cases 
we find a remarkable protein design to transduce a me¬ 
chanical stimulus into a chemical pathway that the cell 
can understand. 

It is much more challenging to understand the mech¬ 


anism of ‘mechanosensitivity of the IP^ kind’, when the 
cell responds to substrate stiffness^^’^^. In order to detect 
a linear response function (elastic constant or viscosity), 
the cell has to apply a force of its own, and two sepa¬ 
rate measurements have to be carried out: of this force 
or stress as it is transmitted to the substrate; and of the 
displacement or strain in the deforming substrate. The 
“stiffness” is a ratio relating these two physical param¬ 
eters. For this kind of measurement, the apparatus re¬ 
quires two points of application: to measure displacement 
a fixed reference point is always needed - to measure force 
one has to ensure a point of reaction as well. We need 
two fingers to squeeze a test object from two sides to de¬ 
termine how stiff it is - in the same way, all engineered 
devices ultimately have two points of action on the test 
sample. A mechanical sensor like the latent transforming 
growth factor-yd complex we study in this paper (Ic-TGF- 
fS) is effectively a single localised molecular protuberance 
outside the cell surface; it only has one point of appli¬ 
cation through which a pulling force is transmitted from 
the cell interior, and no a priori information about the 
reference point in the substrate^^. An attempt to use a 
dynamical regime to overcome the need for a separate 
reference point in space also does not work: it is widely 
known that cells exert a constant force on their envi¬ 
ronment, not a dynamically varying force^^“^^. A recent 
paper by Escude et have formulated a new model, 
using the viscoelastic response of the substrate coupled 
with the stochastic process of rupturing the Ic-TGF-^d by 
a constant pulling force generated by the cell. The result 
of this model is the predicted rate of free TGF-^d release, 
expressed as a compact function of the pulling force and 
the substrate stiffness. 

Gells in the body and in-vitro are normally supported 
biochemically and mechanically by the extracellular ma¬ 
trix (EGM), or its fragments (such as fibronectin) bound 
to artificial substrates. Experiments have observed that 
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the stiffness of a substrate can stimulate differentiation of 
a fibroblast into a myofibroblast, or determine the future 
lineage of stem cells^^: whether they will differentiate to 
bone, muscle, or nerve cells for example. How does the 
cell use proteins to sense stiffness and how does it convert 
a measurement it might be able to make into a chemical 
signal that the cell can understand and respond to? And 
what might a robust system which is designed to answer 
these two questions look like? 

In this work we develop a coupled set of kinetic equa¬ 
tions that describe a dynamical signalling feedback loop 
capable of interpreting the measurement of the substrate 
stiffness and convert it into the morphological changes 
the cell needs in reaction to such a substrate. We fol¬ 
low the earlier work by using the Ic-TGF-/? family as 
a particular biological sensor, shown in^^ to be associ¬ 
ated with the fibroblast to myofibroblast differentiation, 
and in^^’^^ with the development of smooth muscle cells 
from embryonic stem cells, as well as their transforma¬ 
tion following arterial injury or significant deposition of 
mechanically rigid plaque. We combine the result on the 
stiffness-dependent rate of latent complex rupture with 
the analysis of the coupled biological system to derive 
a set of equations to describe the dynamics of the pro¬ 
tein pathway which forms the mechanosensing signalling 
loop. We compare the results and predictions to the 
stiffness-dependent behaviour observed in cell-adhesion 
experiments, as well as in-vivo observations, and find a 
strong correlation of the model behaviour and real sys¬ 
tems. 

In the discussion of the biological system and com¬ 
parison to experimental results we will appeal most sig¬ 
nificantly to the fibroblast system for guidance. This 
is because a substantial volume of research has been 
carried out into the behaviour of this system and be¬ 
cause the model developed in^^, which we incorporate 
into our model, bases its arguments on the fibroblast 
system too. Besides, fibroblasts are the most prevalent 
cells in connective tissue and are vital in tissue repair 
within the body. When a break in the tissue occurs fi¬ 
broblasts differentiate to myofibroblasts which act to pull 
the wound back together. The highly contractile myofi¬ 
broblast form is produced by high expression of a-smooth 
muscle actin (a-SMA) fibres^^, the production of which 
has been shown to be induced by TGF-/?^^’^^. 

The mechanosensing system is summarised diagram- 
mat ically in figure 1. The latent complex of TGF-/? is 
the stretch-dependent mechanosensor; it is representative 
of a larger family of mechanosensor complexes. There 
are three isoforms of the signalling molecule: TGF-/dl, 2 
and 3 with distinct effects in vivo, but their complexes 
are all produced in the same way: tightly but nonco- 
valently bound as a dimer to their latency associated 
peptide (LAP-1,2,3) to form the ‘small latent complex’. 
This forms a larger complex (which we are referring to as 
lc-TGF-/d) by disulphide links between certain domains 
in LAP and the latent TGF-/? binding protein (LTBP), 
which in turn can bind to the EGM to anchor the com- 



FIG. 1. A diagram showing the mechanosensor system, 
from^®. The large latent complex of TGF-/3 binds to ECM 
proteins attached to a substrate, which may be deformable. 
The other end of lc-TGF-/3 is attached to the cytoskeleton 
via an integrin complex that transmits the pulling force. If 
released, the free active TGF-/3 can bind to a set of receptors 
on the cell surface to initiate the signalling loop discussed in 
this paper. 

plex on the substrate^^. The large latent complex as a 
whole is secreted out of the cell, where it may connect 
LTBP to ECM components. Though no unambiguous 
binding partner for the complex has been found in the 
ECM24, there is evidence in favour of strong covalent 
binding, likely to 8-Cys domains on LTBP^^ and it has 
been known to bind to certain ECM proteins such as fib¬ 
rillin and fibronectin^^. We shall assume that this bond 
is strong and not breaking during the processes we study. 

On the other side, the Ic-TGE-/? adheres to integrins on 
the cell surface, which are bound to the actin fibres mak¬ 
ing up the contractile cytoskeleton, allowing the pulling 
force generated by myosin motors to be transmitted to 
the complex, and through it - to the substrate. Upon 
mechanical rupturing of the complex, the active TGE-/? 
is released. Once the Ic-TGE-^d is spent, it can no longer 
function as the adhesion point and may or may not be 
replaced by a new latent complex. 

Once TGE-/d is released, it is free to diffuse. It could 
therefore bind to specific receptors on the cell surface. 
Binding initiates the signalling loop in the cell illustrated 
in figure 3 below. It may stimulate production of new Ic- 
TGF-/3 units to replace those that are spent, and also 
add more a-SMA on the inside of the cell^^. a-SMA is 
produced to form more actin filaments and thus increase 
the pulling force F applied to the latent complex adhe¬ 
sion points. We shall assume a simple linear relationship 
between the concentration of a-SMA and the force F 
applied to Ic-TGE-^d. 


II. MECHANOSENSING REACTION KINETICS 
Rate of latent complex breaking 

To describe the underlying physics of the system we 
start with the rupturing of Ic-TGE-/? by an external 
pulling force F. On the far side this complex is attached 
to a viscoelastic substrate which we model by a simple 
Voigt model with an effective stiffness k, that is directly 
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related to the Young modulus E of the elastic medium via 
the classical continuum-mechanics relation: hz ^ . 

Here ^ is the characteristic length scale of short-distance 
cutoff in the Green’s function of linear elasticity. We 
shall find that our results best match experimental ob¬ 
servations when this length scale is taken as ^ ~ lOOnm, 
which happens to be a characteristic mesh size of the ex¬ 
tracellular matrix of crosslinked filaments. Then a rigid 
glass substrate, with a Young modulus ^ ^ 50 GPa, gives 
the maximum value of our ~ 21 , 000 N/m^. A typ¬ 
ical muscle, organ or connective tissue (as well as syn¬ 
thetic gels) have E ^ 10 — 50kPa, giving the range of 
K ^ 0.004 — 0.02 N/m^. At the other end of rigidity 
scale, the brain tissue is very soft: E ^ 100 Pa, leading 
to the lower bound of /i: ~ 4 • 10“^N/m^. 

Overdamped thermal fluctuations dominate the re¬ 
sponse of this series of mechanical elements. In equi¬ 
librium (i.e. with no force applied) the probability (rate) 
of spontaneous breaking of the latent complex is pro¬ 
portional to exp[—AG/Zc^T] with AG the bonding free 
energy of the large latent complex and ksT the thermal 
energy of the reservoir. When the pulling force is applied 
to the complex, it is also transmitted to the deformable 
substrate and the effective potential holding the latent 
complex together changes. 

The rate of breaking of the latent complex attached to 
a viscoelastic substrate, under a constant pulling force, 
is given by the expression below (which corrects a typo 
in the original work^^): 

;^l/2 A-1/4 -F‘^l2KkBT 

km{E,E)=C-, -^^^-, ( 1 ) 

(A^^ + (iK)l/2^ 1 + AG/feBT 

where the shorthand non-dimensional notations used are: 

^_3Df AG A 2Fwo 

ul \2irkBT) ’ AG ’ ^ SAG ‘ 

Here uq is the characteristic distance the latent complex 
needs to be stretched for rupturing. We could estimate 
the bond free energy AG holding the large latent com¬ 
plex together as for a few hydrogen bonds^^, AG ~ 10 - 
16kcal/mol, and the stretching distance before rupturing 
as the size of an aminoacid residue, uq ^ 0.3 nm. This 
is consistent with descriptions of force induced TGF-/d 
activation by binding of ay^Q to the latent complex^ 
In this case the range of the scaled (non-dimensional) k 
is between 4 • 10^ at the upper bound of rigid glass, 0.02 
for a typical gel or muscle tissue, down to 8 • 10 “^ at the 
lowest bound of = 100 Pa. 

The effective diffusion constant 44, which enters as a 
common factor in the expression for the rate of mechan¬ 
ical breaking km^ is the ratio ksT /reflecting the 
geometric mean of the damping coefficients (loss factors) 
of the substrate and the latent complex (see^^ for de¬ 
tail). We will find it necessary to estimate the value of 
this effective diffusion constant, which is determined by 
the effective damping coefficient of the viscoelastic sub¬ 
strate medium 7 « ktr (that is, the characteristic time 



FIG. 2. Rate of rupture of lc-TGF-/3 (in real units of 
[s“^]) plotted against the non-dimensionless scaled force / = 
Fuo/AG, at various different values of dimensionless elas¬ 
ticity k. The plot highlights the dramatic increase of the 
rate on stiffer substrates, and also the highly non-monotonic 
variation with force. The inset shows that on very soft sub¬ 
strates the rate of rupture is very low and in fact decreases 
further with the force. Note that the rate always diverges at 
fc = 3/2, which is the maximum force that the latent complex 
can stand. 


of stress relaxation in the medium). There is a broad 
range of stress relaxation times found in different ma¬ 
terial systems, but for the sake of simplicity we shall 
take the characteristic Rouse time for a polymer network: 
tr ~ 10 “^s; whether it is a synthetic polymer gel or an 
organic tissue - this may serve as a first estimate (we 
must be aware that such a simple argument is likely to 
fail in the limit of very rigid glass and also in very soft 
gels/tissues with extra-long polymer strands). Then for 
a mid-range substrate stiffness of a typical gel or muscle 
tissue {E = 20kPa), at room temperature we will have 
D « 5* 10“^^m^/s, a sensible value lower than the typical 
constant of thermal diffusion in water. We take the form 
of D to be ksT/hcrji in the numerical analysis below. 

In spite of its apparent complexity, the expression for 
the rate of mechanical breaking of the latent complex 
has only a few straightforward features: it describes the 
latent complex as having the physical bond energy AG 
in equilibrium, which is decreased when a pulling force 
is applied to it, making the basis for the Bell formula for 
the rate^^: when the energy barrier is large, the applied 
force is small, and the substrate stiffness is low - eq. ( 1 ) 
does reduce to: 

km{F, k) « 7 G exp 

However, we cannot use this Bell formula since we need 
to consider high pulling forces and moderate barriers that 
occur in practice. 

On the other hand, eq. ( 1 ) also contains the effect 
of the deformed substrate, which is also subject to the 


AG - Euo 


( 2 ) 
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same force F, transmitted through the latent complex 
while intact. This produces the ‘enzyme effect’ discussed 
in^^: on increasing the pulling force and deforming the 
substrate into a new (deeper) elastic energy miniumum, 
the fluctuations of the point of adhesion become con¬ 
fined; this is expressed by the factor exp[—F^/2 /i:/cbT] 
in eq. (1). The other effect of the substrate viscoelastic 
response is to alter the non-exponential prefactor of the 
rate expression by adjusting the statistical weight of the 
intact latent complex in the full ensemble. The depen¬ 
dence of the breaking rate, on the applied force for 
several values of substrate stiffness is illustrated in figure 
2 . 

At sufficiently high substrate stiffness, expressed by 
the non-dimensional scaled parameter k given below the 
eq. (1), there is a maximum in the rate of breaking of 
the latent complex. So, initially, the positive feedback 
is expected - an increase in the pulling force increases 
the rate of latent complex breaking, leading to a further 
increase in force, etc. until the cell finds the equilibrium. 
However, at very low values of k the maximum is at zero 
applied force. That is, on such soft substrates a pulling 
force actually reduces the rate of latent complex breaking 
causing the negative feedback. The outcome depends on 
the details of interactions between the key parameters in 
this feedback loop. 


Protein interactions 

Upon rupturing of one latent complex, an active TGF- 
P unit is released and becomes free to diffuse towards a 
specific receptor type on the cell membrane, and bind to 
it. The binding of TGF-^d stimulates production of a- 
SMA^^ by a signalling pathway that we do not go into 
details of, except to know the ratio of how many a-SMA 
monomers are produced to each TGF-/3 bound. Given 
the linear relationship between a-SMA and force F ap¬ 
plied to Ic-TGF-yd, this stimulated production results in 
an increased force on the latent complex. The increase 
in applied force changes the rate of rupture of the com¬ 
plex, with the feedback nature depending on the regime 
of km{F,K). 

We now seek to mathematically model the system of 
protein interactions described above. The three key pro¬ 
teins are a-SMA (labelled a), TGF-^d (labelled and 
Ic-TGF-yd (labelled A), see figure 3. To be specific with 
dimensionality, we consider the evolution of the number 
per cell of each of these proteins, a{t)^ P{t) and A(t). The 
number per cell can be re-dimensionalised into a concen¬ 
tration by considering the typical mass of the protein 
in question and the volume of the cell - e.g. a typical 
fibroblast. 

Within our model there are certain assumptions about 
the state of proteins that we count in our variables. A(t) 
represents only those lc-TGF-/3 that are attached to both 
the cell and the substrate and have not ruptured. P{t) 
includes only TGF-^d molecules that are released from 
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FIG. 3. A diagram showing how the proteins in our 
mechanosensing model interact. The three key markers of 
the signalling loop are: A (the concentration of lc-TGF-/3), [3 
(the concentration of active TGF-/3) and a (the concentration 
of SMA). The rate of spontaneous breaking of A, /cm(o, F), is 
a non-linear function of the pulling force F applied by the 
cell, eq. (1). 

their latent complex and are close to the cell surface, 
but not yet bound to a receptor. We can describe the 
evolution of these three variables in the following set of 
differential equations: 

^ =-km{a,K)X +Axkp/i (3) 

^ = +km{a, k)X - kpj3 (4) 

— = X-Aakjip. (5) 

Here the coefficients Ax and are, respectively, 
the number of fresh lc-TGF-/3 complexes and a-SMA 
monomers produced for every instance of TGF-/3 binding 
to a receptor. Implicit within these two quantities are all 
the underlying protein pathways whose starting point is 
a TGF-/3 binding to its receptor, and whose end point 
is the production of a new latent complex and a-SMA 
subunit, respectively, kj^ is the rate of binding of TGF-/d 
to the cell receptors. The key element of this model is 
the rate of the rupture of Ic-TGF-^d, given in eq. (1) and 
discussed above. 

We have to assume a linear relationship between a{t) 
and the pulling force F to be of the form: 

F = Cfmnma, ( 6 ) 

where is the fraction of a-SMA monomers polymerised 
into actin filaments that are involved in the application 
of force on A (this is expected to be a small number). 
is the force exerted by a myosin motor, which is a known 
parameter for each individual motor (^ 2pN), but may 
be a function of ATP availability in the vicinity of the 
adhesion complex; rim is the number of myosin motors 
that act along each filament on average, a parameter that 
the cell can adjust within a wide range. 

The dissipative sinks 

We must also consider certain dissipative terms that 
affect all three of our dynamic variables, so that the sys- 
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tern can find a homeostatic state; three dissipative terms 
are relevant, one for each of the three proteins. 

The A sink is the saturation of integrin binding sites. 
We expect there to be only a certain number of sites over 
the surface of cell in contact with the substrate, at which 
the latent complex can attach to both the cytoskeleton 
and the substrate. This number may be large, but it 
is nevertheless limited, and this restricts the maximum 
value that A can take. We consider a probability of a Ic- 
TGF-yd attaching to the substrate and cytoskeleton which 
has the form: 

P{X,N) = Po^^, (7) 

where N is the total number of available binding sites on 
the cell surface and Pq is the probability of a latent com¬ 
plex attaching to the substrate and cytoskeleton when 
N = oo (i.e. there is no limitation on A). This proba¬ 
bility should be multiplied by the rate of production of 
lc-TGF-;d {Axkj^p) to give a total rate of production of 
A, to replace the second term in eq. (3): 


dj\ A 1 ^ ^ n 

M=AxhPo^;^P. (8) 

Given the context of the problem, we expect the free 
TGF-yd molecules to not only diffuse towards the cell re¬ 
ceptors but also could drift away from the cell entirely. 
Therefore, a (3 sink due to diffusion away from the cell 
surface should be included. We consider this diffusion 
loss by treating TFG-^d as a free Brownian particle dif¬ 
fusing away from a sphere (representing the cell surface). 
By considering the Gaussian probability for a molecule 
to be at a position r from the surface after time t, then 
setting r = 0 we find the concentration near the surface 
to decrease as P{r = 0) = ^ where £ is some 

normalisation constant with units of length. 

Our variable P{t) is then considered to be the proba¬ 
bility of finding the TGF-^d particle at r = 0. Initially 
we expect TGF-^d to be present only at the cell surface, 
where it was released from the latent compolex. To nor¬ 
malise to unity at t = 0 we need to identify the radial 
lengthscale on which TGF-^d exists. Because the proba¬ 
bility tends to a delta function as t tends to 0, we assume 
that the radial lengthscale for normalisation is the diam¬ 
eter of TGF-;d molecule: £ = 2i?^. 

Gonsidering the probability as a volume fraction, we 
can write P(0, t = 0) = ydV^/Vtot where is the volume 
of one TGF-P subunit and Vtot is the total volume oc¬ 
cupied by all /3 at t = 0: a spherical shell at the radius 
of the cell (i^ceii) of thickness 2i?/3. Therefore we want 
the contribution to the differential equation (4) to yield: 
(3{t) = 2RpVtot/{VpV‘^7TDt), when no other terms are 
present. Differentiating gives a rate of diffusive loss: 


dl3 _ nDV! 

dt 


(9) 


For our calculations it will be important to estimate the 
coefficient in this relation, which involves several known 


length scales; the values for D, Rcew and R^ are given in 
table below. 

In the cell we always find that a-SMA decays over 
time, being used in various other processes. The rate of 
this spontaneous decay has been measured experimen¬ 
tally and also listed in the table below; we denote it as 


= -kaa. ( 10 ) 

Gombining all of the above together we arrive at the 
final set of three rate equations, which are now have bal¬ 
anced sources and sinks for the evolution of all three of 
our markers: 


dA 

dt 

dp 

dt 

da 

dt 


N - X 

^mipi “h jAxPohp ^ jd 

( 11 ) 

ttDV! , 


+ km{a,K)\ kpp 2 t/2 

^-^heta ^tot 

( 12 ) 

“h Afj^k P kf^a 

(13) 


Fixing parameter values 

The non-linear kinetic model of the signaling feed¬ 
back loop has many parameters. Where available, values 
of experimentally derived parameters were always used in 
preference. The table below gives the values of parame¬ 
ters that indeed have been determined experimentally. 


Experimentally determined values 


Parameter 

Value 

Origin 

-Pcell 

5 fim 

31 

k(x 

6- 

30 

fm 

2.3N 

32 

Tim 

80 

33 

N 

1000 

34 

D 

7.7- 

35,36 

Rp 

1.66 nm 

37 


The value of k^ was obtained from the line in figure 5 of^^ 
representing decay of a-SMA without TGF-^d in a cell an¬ 
chored to a gel. By assuming exponential decay between 
day 1 and day 5, a value was approximated which fits this 
curve. The value of D is the diffusion constant for Green 
Fluorescent Protein (GFP), which is comparable in mass 
to TGF-;d, and so it is expected that the values of both 
diffusion parameters are very close. No value of Rp could 
not be found in the literature so it was approximated by 
assuming that GFP (a cylindrical protein) and TGF-/3 
have similar volumes and that TGF-/3 is spherical. The 
diameter and length of GFP are 2.4nm and 4.2nm re¬ 
spectively. The values of D, AG and uq for the large 
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latent complex mechanosensor itself have been discussed 
earlier in the text. 

Other parameters were not experimentally defined, so 
values that seemed physically reasonable were chosen 
and the dependence of system behaviour on these val¬ 
ues checked. 


Model assumptions 


Parameter 

Value 

PqAx 

1.7 

Act 

3 

c 

I • 10“^ 


The value for PqAx was selected by assuming each 
TGF-/3 binding to a receptor stimulated production of 
approximately 1 or 2 new Ic-TCF-^d. We expect all (or 
almost all) the Ic-TGF-yd produced to be used in this 
mechanosensing system. We also assume that Pq is close 
to 1. A value of 1.7 was selected to approximately reflect 
these assumptions. 

Aa was selected to roughly reflect the high concentra¬ 
tion of a-SMA in the cell. In order for the cell to react 
to a stimulus it will need to produce a significant amount 
of a-SMA monomers. 

The value of C, the fraction of newly produced a-SMA 
that ends up polymerised into F-actin and contribute to 
the pulling force on the adhesion complex, was chosen by 
considering the ratio of the volume of an a-SMA filament 
monomer (modelled as a sphere) to that of an a-SMA fil¬ 
ament subunit (modelled as a cylinder). We define the 
length of a subunit as the step size of myosin which is 
8 nm^^. The diameter of actin filaments is about 8 nm^^. 
The average diameter of an a-SMA monomer is about 
4.8 nm^^. This gives a ratio of 0.14. We then multiply 
this ratio by the fraction of a-SMA being used in this 
process to total a-SMA in the cell. This was estimated 
roughly to be 1/1500, but whatever this value - it is 
reassuring that this number is small since one does ex¬ 
pect that only a small fraction of new a-SMA monomers 
would contribute to the increase of the pulling force on 
the point of adhesion. 

The most difficult parameter for us to estimate is the 
rate of binding of free TGF-P to its receptors on the outer 
cell surface, kj^. It is not experimentally determined and 
is difficult to define through arguments of plausibility. 
We fixed its value by requiring that the transition ob¬ 
served in steady-state values of A, and a (discussed 
in the following sections) corresponds to experimental 
results, so in a way is the only true fitting param¬ 
eter of our model. In culture the fibroblast transition 
occurs at a Young modulus value E ^ 11 kPa^^. Us¬ 
ing n = (as above, with ^ 100 nm) we find 

the corresponding dimensionless (scaled) elastic modu¬ 
lus: ^transition ~ 5.9 • 10“^. To rcproducc this value in 
our results, the rate needs to be in the range: 


Constrained parameters 


Parameter 

Value 

ki3 

5- 


Method of solving. 

All studies of the differential equations were done in 
Matlab using the ‘ode23s’ ordinary differential equations 
numerical solver package provided. Initial conditions 
were always set randomly, unless otherwise stated. For a 
given set of initial conditions the system always relaxed 
to the same steady state values. For k > 1 the solver 
would often take a very long time to produce a result, 
depending on initial conditions, because of the extreme 
stiffness of the system. Combinations of initial conditions 
and k that took a long time to solve were avoided for the 
majority of the calculation. However, the behaviour was 
checked with one or two ‘awkward’ combinations to en¬ 
sure that it agreed with other results and trends given in 
this paper. Note that to ensure A reached steady state, 
a very long time span had to be used. 


III. EVOLUTION OF VARIABLES 

The variables a, f3 and A evolve in time from various 
initial conditions to the final steady state. Evolution 
maps of the variables are shown in figure 4 for low k 
and a reasonably high k. The value for which we con¬ 
sider k to be ‘high’ or ‘low’ will be discussed later. It is 
clear that the dynamical system is robust and has a sin¬ 
gle stable attractor for a wide range of parameters (we 
have tested a large number of parameter values, not just 
the ones given in the tables above). The position of this 
fixed point attractor is dependent on k, that is, a single 
fixed point for each value of k, and therefore provides a 
vehicle for mechanosensing of the IP^ kind. 

Figure 5 shows the dependence of the time taken for 
f3 to reach steady state (relaxation time) on k. On the 
whole, on soft substrates the relaxation time varies from 
about 10^‘^ s to a little over 10^ s (8 to 30 hours). Ex¬ 
periments have shown that cells react to TGE-^d on time 
scales as long as 48 hours (upper bound^^’^^) - agreeing 
well with our results. It is interesting to observe how the 
time taken to reach the equilibrium steady state at high 
k increases substantially. At maximum possible values of 
substrate stiffness {k > 100) the average value is about 
10^-^ s (or 3-4 days). It is possible that this high val¬ 
ues of k are not practically achievable at all since rigid 
materials with a Young modulus E ^ I-IOGPa have, in 
fact, a much smaller internal mesh size However, the 
overall trend of increasing the time to reach the homeo¬ 
static equilibrium with the substrate stiffness is a definite 
prediction of our theory. 
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FIG. 4. Graphs showing the evolution of a, (5^ and A on 
logarithmic scales for: (a) low value of substrate stiffness k = 
10“®, and (b) high substrate stiffness k — 0.1, from different 
initial conditions. The final steady state reached for (a): A = 
412, /3 < 1 and a — 1.6, and for (b): A < 1, /3 = 108 and 
a = 2718. 


IV. THE BEHAVIOUR IN THE STEADY STATE 

The steady state dependence of each variable on k is of 
the highest importance in comparing our model to real 
biological systems and seeing how well it can describe re¬ 
ality. In a real cell, we would expect the concentration of 
TGF-/3 to determine the response of the cell. Different 
concentrations would result in different gene expression 
or different enzyme activation and overall a different be¬ 
haviour by the cell. We use the fibroblast system as the 
basis for discussion in this section. 

Figure 6 shows how the concentration of each variable 
behaves in the steady state. There are three distinct 
regions: the low and high k regions where all variables are 
unchanging, and a middle k region which is the transition 
region. 

We first consider some peculiarities of the behaviour 
of A. In the high k region we see the value of A tends to 
very close to zero, suggesting that as soon as lc-TGF-/3 



FIG. 5. A graph of the time for each simulation to reach the 
equilibrium steady state, defined as when /3(t) reaches its equi¬ 
librium value within a tolerance of 1 per cell, plotted against 
the effective substrate stiffness k. Each point represents an 
independent simulation with randomly selected initial condi¬ 
tions, varied by four orders of magnitude. The arrow marks 
the experimental position of the myofibroblast transition and 
the dashed lines mark the ‘crossover region’ of Fig. 6. 



FIG. 6. A graph showing variation of concentration of steady- 
state values of a, /3 and A (in units of particle per cell) against 
thescaled substrate stiffness k (expressed in decimal logarith¬ 
mic units). Initial conditions were random. The dashed lines 
show the approximate location of the transition region and 
the arrow marks the position of myofibroblast transition. 


binds to the integrin and substrate - it ruptures, so that 
on average there is none (or very little) present. This is a 
strange result and is certainly not intuitive. The reason 
for this is that in the high k region focal adhesions form 
(illustrated in figure 7) where many lc-TGF-;5 cooperate 
together in a combined structure. This is not accounted 
for in our model because we assumed that A represented 
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FIG. 7. A diagram of a contracted myofibroblast and less con¬ 
tracted fibroblast. Adapted from a discussion in^^. The large 
increse in ct-SMA expression is evident (as is widely reported 
increase in active TGF-/3 on stiff substrates) - but our conclu¬ 
sion about A (linked to the other two markers) suggests that 
there are only large focal adhesions, and no/few individual 
lc-TGF-/3 on the cell surface on stiff substrates. 



only lc-TGF-;5 attached individually to the substrate and 
integrin. Therefore we can deduce that only a small frac¬ 
tion of total lc-TGF-/3 exist as individual elements in the 
high k limit - a deduction supported by experiments^ 
showing that focal adhesion size increases on stiffer sub¬ 
strates. The true number of lc-TGF-/3 may be much 
higher than that implied by the value of A. 

Examining the behaviour of a and (3 in steady state, 
we see that the transition occurs in the region 10“^ < 
k < 10“^-^. The corresponding Young’s modulus where 
the fibroblast transition occurs lies in between these two 
extremes, as required by our fitting parameter. 

Over this transition p increases from 6 cell” ^ to a final 
value of 108 cell” These values correspond to real con¬ 
centrations of about 0.48ng/ml and 8.6ng/ml, respec¬ 
tively - taking TGF-/5 molecular mass as 25kDa^^ and 
assuming the fibroblast is a spherical cell with a radius 
of 5 /im (a big approximation given the myofibroblast is 
strongly contracted). How does this compare to real sys¬ 
tems? TGF-;5 concentrations as low as O.lng/ml were 
found to be sufficient to induce a transition in endothe¬ 
lial cells to myofibroblasts with a maximal response at 
Ing/ml^^. Another article^^ observed stimulation of the 
fibroblast to myofibroblast transition for concentrations 
of TGF-/3 between O.lng/ml and 15ng/ml. 


FIG. 8. Plot of km{f,k) (in real units of [s”^]) for three dif¬ 
ferent values of scaled substrate stiffness k. The red cross 
represents the value of the dimensionless scaled force /eq in 
cell equilibrium, corresponding to the steady state value of 
a. In the lowest plot (for stiff substrate), the red cross corre¬ 
sponds to a value of km = 0.01 s”^. 

are three distinct types of behaviour. At very low k we 
see the steady state force is beyond the maximum in km 
which is at / ~ 0 (so the cell is not predicted to reduce 
its pulling to zero, as might have been assumed). At 
high k on rigid substrates, the steady state force occurs 
before the maximum, so in spite of its restructuring and 
formation of more stress fibers - the coupled signalling 
loop stops the cell short of reaching the maximum rate 
of latent complex breaking. At intermediate values, in 
the ‘transition region’ of k, the steady state force rests 
somewhere after the peak of km^ 

We find that the transition in A from a high steady 
state concentration at low ^ to a low concentration at 
high k occurs between the low and intermediate cases 
discussed in this section. No particular effect in A is 
observed when the steady state value of a moves past the 
maximum in km- As far as the system has been tested, 
it seems that the maximum in km has very little part to 
play in the dynamics of our coupled system variables. 


The breaking rate km{f^k) and applied force in equilibrium 

Finally, we examine the relation of force and rate of 
rupture of Ic-TGF-^^. Qualitatively it would seem rea¬ 
sonable to expect the steady state to correspond to the 
maximum in km (Fig. 2) where the positive feedback 
changes to negative. Surprisingly we do not find this 
intuitive relation between system variables and km- 

Figure 8 shows the relation between the steady state 
value of a (directly proportional to force) and the rate 
of rupture of Ic-TGF-/?. As can be seen, system vari¬ 
ables do not relax on the maximum in km- Instead there 


Using Hm to adjust the transition region. 

The system proposed is clearly not sensitive to the sub¬ 
strate stiffness in the low or high k regions. However, it 
has been observed that cells respond to elasticities out 
of the range our model would suggest is possible^^. We 
found that changing Um does not affect the steady state 
values of any variables away from the transition region. 
But it does affect the position of the transition region ver¬ 
sus ^ - as seen in figure 9. If we allow the cell the flexibil¬ 
ity to vary Um as well as the applied force, then Um could 
be adjusted so that the transition region scans over many 
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FIG. 9. The effect of varying values of rim (the ‘number’ 
of active myosin motors on the actin fiber) on the steady- 
state values of our three marker variables. The main effect 
is in adjusting the position of the transition region of high 
sensitivity to stiffness. The values are printed on the plot: 
rim = 50 (red), rim = 80 (blue), and rim = 300 (black). 


orders of magnitude of k. Myosin motor dynamics can be 
regulated by de/phosphorylation in-vivo"^^. Therefore, as 
long as the cell has a method to record rim and one of 
the variables, then it could dynamically vary rim to find 
the elasticity of its substrate. 


V. CONCLUSION 

This work was a natural continuation of^^ on the rate 
of rupture of lc-TGF-;5. Incorporating their model into a 
full dynamically coupled system allowed us to create a set 
of differential equations with the mechanical sensing sys¬ 
tem encoded in this rate. The behaviour of the equation 
is entirely dependent on the applied force and the elas¬ 
ticity of the substrate, and we have checked by start¬ 
ing from a large variety of random initial conditions that 
the equilibrium (steady state, homeostasis) of the cell is 
unique and specific on each substrate. We found that 
the system always returns to the steady state indepen¬ 
dently of initial conditions (a single stable attractor) and 
there appears to be a certain characteristic time scale of 
relaxation towards this equilibrium, regardless what the 
cell parameters have started with, of the order of 10^‘^s 
(i.e. 8-9 hours), which is a broad agreement with in-vitro 
observations of cell deposited on various substrates. 

The choice of model parameters may appear the most 
controversial aspect of our work (since there are so many 


of them and some are little-known), yet we have taken 
great care to fix most of these parameters at specific val¬ 
ues found in experiment. It is clear that for any partic¬ 
ular cell system (osteoblast, fibroblast, smooth muscle, 
glial - or indeed stem or various progenitors) the values 
of these parameters would be subtly different, and the 
cell response range vary accordingly. But the range of 
this variation is not great, and we are assured that the 
qualitative nature of the predicted response will remain 
the same. 

To summarise the predictions of our theory: 

1. We study a response of a ‘model cell’ when it is 
placed on a specific substrate, with given elastic modulus 
and loss factor, or in other words - stress relaxation time. 
The principal sensor is the large latent complex of TGF- 
/3, which we use as a representative model. Following^^, 
we found the compact expression for the rate of latent 
complex rupture, km{F, k) in eq. (1), which is deter¬ 
mined by the parameters of the substrate and the steady 
pulling force from the cell. This breaking rate has a com¬ 
plex non-monotonic behavior with changing parameters, 
and it is central to the whole process of mechanosensitiv- 
ity of the IF^ kind. 

2. On release of the signal molecule (active TGF-/3) 
after the irreversible breaking of the latent complex, it 
could bind to a specific receptor on the cell to initiate the 
feedback loop - or it could drift away from the cell. We 
have not yet followed an obvious conjecture, that these 
‘lost’ TGF-/3 signal molecules may provide an element of 
quorum sensing in a confluent tissue or cell colony. The 
concentration of free TGF-/3 released from the broken 
latent complexes (our variable /3) has specific values in 
homeostatic equilibrium. 

3. On placing the cell from any previous environment 
onto the given substrate (which we modelled by scan¬ 
ning a broad range of initial conditions), the cell remod¬ 
els itself until it achieves the ‘equilibrium’ morphology 
specific for this substrate. The most obvious and exper¬ 
imentally tested signature of this is the concentration (3 
in the vicinity of the cell, which is low/vanishing on very 
soft substrates, and high constant on mechanically rigid 
substrates (which is what has been registered in many ex¬ 
periments on very different cells and environments). The 
transition region in between these two limits is where the 
cell sensitively responds to small changes in stiffness by 
altering the equilibrium values of its parameters. 

4. The adjustment of cell morphology to its substrate 
stiffness happens via the expression of additional a-SMA 
(stimulated by the TGF-/5 binding), which is then incor¬ 
porated into additional F-actin and stress fibers - and 
provides more strong pulling force F{a) on the latent 
complexes in points of adhesion. Again, in close agree¬ 
ment with observations, we find that actin production is 
low on soft substrates and high on rigid substrates, with 
a transition region in between. The fibrosis in tissues is 
therefore linked to the stiffness of the environment the 
fibroblasts or smooth muscle cells are in contact with. 

5. Our final kinetic variable, the concentration of 
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intact large TGF-;5 latent complexes (A) shows an inter¬ 
esting an unexpected response to stiffness change. On 
soft substrates, with low a (so the cell does not pull too 
hard on its adhesion points) and low (3 (reffecting a low 
probability of Ic-TGF-^d breaking) - there are a lot of 
adhesion points around the cell: the equilibrium level of 
A is high. However, on rigid substrates when the cell 
pulls with a large force, we find A dropping to a very low 
value. Note that since the level of [3 is high in this region, 
the breaking of latent complexes occurs very frequently - 
almost immediately after they have been replenished by 
the cell, so there are almost no active adhesion points in 
operation around the cell surface. This may sound con¬ 
tradictory, but we interpret it as the commonly observed 
state when just a few large focal adhesion complexes hold 
the cell on a rigid substrate - and we have now discov¬ 
ered that the cell surface outside of these few focal adhe¬ 
sions is free from the ‘singular’ adhesion points. This is 
in contrast with the same cell on a soft substrate when 
it has a homogeneous concentration of singular adhesion 
points (high-A state). Our theory does not describe the 
collective effects that take place in large focal adhesions: 
certainly both the pulling force and the latent complex 
rupturing are very different there. 

We have tried to make sure the stiffness transition re¬ 
gion where steady state concentrations of a and [3 rise by 
many factors, and A falls significantly, corresponds to ob¬ 
servations for fibroblasts. Since most of the parameters 
we use are independently fixed, the critical role is played 
by the binding rate of TGF-^d to its receptor which 
we assume is not known. The value of kj^ ~ 5 • 
was required to bring this transition region to the Young 
modulus E ^ 10-15 kPa required to observe fibroblast to 
myofibroblast transition. However, we have shown that 
the cells have several mechanisms to adjust this transi¬ 
tion region even for exactly the same Ic-TGF-yd adhesion 
parameters: for instance by varying the number and the 
efficiency of myosin motors (the parameter Fig. 9), 
or by channeling a different amount of newly produced 
a-SMA to the adhesion pulling filaments (the parame¬ 
ters C and ka)- Therefore we would expect that glial or 
neuron cells would ‘move’ their stiffness transition to a 
much lower Young modulus range, while osteoblasts may 
benefit from their region of primary sensitivity to lie at 
a higher Young modulus as appropriate to their target 
environment. Perhaps this choice is made at the stem 
level and is at the core of subsequent differentiation: by 
dynamic variation of rim by de/phosphorylation the cell 
can ‘scan’ the transition region over many orders of mag¬ 
nitude of elasticity. 

We found that the system had no particular depen¬ 
dence on the prominent and characteristic peak of the 
lc-TGF-/d rupture rate km{F, /^), defying our intuitive ex¬ 
pectation for how the system might respond. It seems 
that the general pattern of behaviour in homeostatic 
equilibrium is independent of many parameters - at least 
to the extent that it was tested. Further work may seek 
to answer some of many questions posed by this paper. 


Does it matter that the peak in km seems to not have the 
importance we expected? How could we include the ef¬ 
fect of focal adhesions on Ic-TGF-^d and would it change 
the behaviour of the model dramatically? 
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